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Recently the authors developed a scattering approach that allows for a complete description of 
the steady-state physics of quantum-impurities in and out of equilibrium. Quantum impurities are 
described using scattering eigenstates defined ab initio on the open, infinite line with asymptotic 
boundary conditions imposed by the leads. The scattering states on the open line are constructed 
for integrable quantum-impurity models by means of a significant generalization of the Bethe-Ansatz 
which we call the Scattering Bethe-Ansatz (SB A). The purpose of the paper is to present in detail 
the scattering approach to quantum-impurity models and the SBA and show that they reproduce 
well-known thermodynamic results for several widely studied models: the Resonance Level model, 
Interacting Resonance Level model and the Kondo model. Though the SBA is more complex than 
the traditional Thermodynamic Bethe Ansatz (TBA) when applied to thermodynamical questions, 
the scattering approach (SBA) allows access to an array of new questions that cannot be addressed 
otherwise, ranging from scattering of electrons off magnetic impurities to nonequilibrium dynamics. 



I. INTRODUCTION 

Recent advances in nanotechnology have allowed ex- 
tensive experimental study of quantum impurity systems 
out of equilibrium in controlled, tunable settingsA The 
impurities are typically realized as quantum dots, tiny 
islands of two-dimensional electron gas attached to leads 
via tunnel junctions. The number of electrons on the 
dot can be controlled using a gate voltage since the hop- 
ping of electrons is impeded by a large charging energy 
U. When there is an odd- number of electrons on the 
dot the upper-most energy level contains only a single, 
unpaired electron, which behaves effectively as an Ander- 
son or Kondo impurity coupled to the two (or more) leads 
playing the role of electron baths. Applying a potential 
difference between the leads results in a current flowing 
across the dot. A wealth of new experimental data has 
been collected in recent years on quantum-impurities out 
of equilibrium in this setting including current vs voltage 
curves and nonequilibrium density of states (DOS) on 
the quantum dots^. Nonetheless, a comprehensive the- 
oretical understanding of the physics of these models is 
lacking. 

Quantum impurity systems are also the simplest ex- 
amples of strongly correlated electron systems, wherein 
interactions between electrons are strong enough to re- 
sult in new collective behaviors which require a new set 
of degrees of freedom for their description- the Kondo ef- 
fect being a canonical example^. The strongly-correlated 
behavior is typically characterized by a low energy scale 
such as the Kondo temperature below which strong cor- 
relation physics dominates and perturbative descriptions 
break down. One of the most fascinating new frontiers 
in strongly-correlated systems is the study of such sys- 
tems in out-of-equilibrium situations. Quantum impu- 
rities are an ideal experimental and theoretical setting 
for exploring the interplay between nonequilibrium- and 
strongly-correlated dynamics due to the relative simplic- 
ity of these models and the wealth of experimental data 
available. 



New theoretical questions arise in this context. Do suf- 
ficiently large voltages suppress strong-correlations and 
thus kill the Kondo effect? Do new scales, such as the 
decoherence scale, arise? Does voltage effectively behave 
as a temperature? How should one handle intrinsically 
non equilibrium phenomena such as nonequilibrium par- 
ticle and energy currents or entropy production. What is 
the effect of strong correlation the entropy production? 

Currently, the most commonly used technique to treat 
quantum-impurities out of equilibrium is Keldysh per- 
turbation theory The perturbative methods, however, 
are applicable only in the high voltage regime and break 
down precisely where strong correlations become impor- 
tant. As such, they are unable to answer the interesting 
questions proposed above. A variety of non-perturbative 
techniques have been developed in order to capture the 
strong correlation physics of quantum impurity models, 
mainly in the context equilibrium physics. These in- 
clude renormalization group methods, techniques from 
bosonization and conformal field theory, and exact so- 
lutions using the Bethe-Ansatz'^. Most of these methods 
arc no longer applicable when the influence of nonequilib- 
rium dynamics is comparable to the strong correlations 
in the problem. This highlights the need for new the- 
oretical approaches that can probe the interesting non- 
perturbative regimes^iiiSiiSiiiS. 

Recently we have introduced such a non-perturbative 
framework that allows the description of a steady state 
out-of-equilibrium quantum impurity system in terms of 
a time-independent scattering formulationii. A steady 
state ensues^ when the system is open. Open systems 
must be defined directly on the infinite line to allow an 
in- flow and out-flow of electrons and energy from the sys- 
tem. The infinite volume limit, which needs to be taken 
ab initio, provides a dissipation mechanism. Under these 
circumstances the non-equilibrium steady state can be 
described by a scattering eigenstate of the full hamilto- 
nian, an eigenstate defined on the infinite line with its 
asymptotic behavior specified at the incoming infinity. 
In most cases the asymptotic boundary conditions are 
determined by the electron leads^i. 
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We have subsequently also introduced a method, the 
Scattering Bethe Ansatz (SB A), to construct those scat- 
tering eigenstates on the infinite line for the Kondo model 
and other integrable impurity models. The traditional 
Bethe Ansatz, on the other hand, which has been exten- 
sively applied to these models, is defined with periodic 
boundary condition with periodicity L (with L subse- 
quently sent to infinity) . This approach is appropriate to 
closed systems and allows an efficient calculation of the 
thermodynamic properties of the systems. However, it 
does not give access to their scattering properties, nor to 
the non-equilibrium physics. 

The scattering approach can also be applied under 
equilibrium conditions, when all baths are held at the 
same chemical potential, or in the case when only one 
lead is present. The purpose of this paper is to study 
the Scattering Approach under these simpler circum- 
stances and confront it with the conventional approach 
which can also be applied here. We will show that 
the SBA reproduces known thermodynamical results for 
quantum-impurity models. Nonetheless, as mentioned 
above, the algebraic Bethe-Ansatz and its finite tem- 
perature counterpart the Thermodynamic Bethe Ansatz, 
prove technically easier when calculating thermodynam- 
ics of quantum-impurity models. The real advantage 
of the SBA is that it allows us to harness the power 
of integrability to explore new questions about electron 
S-matrices and T-matrices important for understand- 
ing quantum-mechanical coherence and dephasing due to 
magnetic impuritie o^"^'^^ . And, in a context which we will 
not further explore in this paper, SBA allows us to un- 
derstand nonequilibrium steady-states in these models. 

The paper is organized as follows. We start with a for- 
mal introduction to the scattering approach to quantum- 
impurity models and the scattering Bethe-Ansatz (SBA). 
Subsequently, we demonstrate our ideas on the Resonant 
Level model where the physics is particularly simple since 
the Hamiltonian is quadratic. Finally, we use the SBA 
to construct scattering states to reproduce well-known 
T = equilibrium results for Interacting Resonant Level 
and Kondo models. We conclude the paper with some 
conjectures about the Kondo model that significantly 
simplify the calculation of some impurity properties. 



II. THE SCATTERING FORMALISM 

The basic idea underlying the scattering formalism is 
the observation that a quantum impurity can be viewed 
as a localized dynamical scatterer off which electrons 
from the attached leads or host metal scatter. The scat- 
tering changes the internal state of both the impurity and 
host electrons and leads to the generation of strong corre- 
lations. The standard procedure for treating such prob- 
lems is to set-up an initial state with wavepackets that 
represent the incoming particles in the far past and evolve 
this initial state for a very long time with a time-evolution 
operator U{t,to) = Texp(— z H{t')dt') of the appro- 



priate interacting field-theory; H is the Hamiltonian that 
describes the particles and dynamical scatterer - in this 
case the quantum impurity. As the impurity is local, the 
interaction switches off far away from the impurity and 
we can define 'in' or 'out' states by specifying the asymp- 
totic behavior in the far past or the far future. Namely, 
the 'in' , respectively, 'out' states are eigenstates of the 
total Hamiltonian, satisfy the boundary conditions that 
they tend to plane waves representing free incoming par- 
ticles in the t — > — cx) and t — *■ cx) limits respectively. The 
cross-section for a particular process is then obtained by 
calculating the overlap of the "in" state with an appro- 
priate "out" state. A recent apphcation of these ideas 
to quantum impurities is given i n^^i^^ . While scatter- 
ing states are designed to allow access to the scattering 
properties of the system, they also allow calculation of 
the thermodynamic properties. 

The Hamiltonian for a quantum impurity attached to 
a bath of free electrons is of the form: 

with efe the full three-dimensional dispersion of the elec- 
trons and a denoting the internal degrees of freedom of 
the electrons. The term i/jnt describes the impurity and 
its interaction with the bath of electron. Examples in- 
clude the Kondo interaction, 

ak a'k' 

with S a localized spin representing the impurity, and 
the resonance level model (RLM), 

Hint ^iY.^'4'l^ + h.c.) + tddU, 

k 

describing a local level at energy e^j which hybridizes with 
the bath electrons. 

Standard manipulations^'' allow us to rewrite the 
Hamiltonian as chiral 1-d field theories. Since only the 
combination r "0^ - enters into the interaction we can 

^ak 'f'ak 

rewrite the theory in terms of the field 

= J d^k Siek - 
as {D denotes the bandwidth, namely the cut-off) 

.D 

Ho= dee V'LV'cc 

J-D 

while the interaction terms take the form, 
J / deipl^ i^)aa' J de'ipa-.e'-S OT t{J de ipld+ h.c.) + edd^ d 
for the Kondo or the IRLM Model respectively. Finally 
introducing a chiral fermion field 
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the kinetic term becomes 

Hq = -i / dx i^iix)dxipa{x) 

while the field enters locally into Hint, in the form "01(0)- 
As we are interested in the physics on energy scales much 
smaller than the cut-off we consider only universal 
results obtained in the limit _D — > cx). 

Open vs Closed Boundary Conditions: The 
scattering approach to quantum impurity problems, by 
its very nature, is defined in infinite systems, without 
boundaries. Physically, this is equivalent to requiring 
that once incoming electrons scatter off the impurity they 
do not return and scatter off the impurity again. We re- 
fer to infinite size systems with no boundaries as "open 
systems". The infinite size of the electron bath assures 
that the host metal or lead is a good thermal bath. Real 
life systems are not infinite but possess boundaries; our 
treatment is valid as long as the return time for the elec- 
trons is much smaller than the system size. The infi- 
nite size of the system implies that scattering states are 
no longer normalizable and in particular, the Feynman- 
Hellman theorem no longer holds^^. This will be impor- 
tant in understanding the results of later sections when 
we construct eigenstates for the IRLM and Kondo mod- 
els. 

In this open system framework the nature of the in- 
coming particles that scatter off the impurity is specified 
by asymptotic boundary conditions. The incoming par- 
ticles, far from the impurity, are eigenstates of the free- 
electrons Hamiltonian Hq and any eigenstate of Hq is a 
possible boundary condition describing what the incom- 
ing particles look like. Two different boundary conditions 
are of primary interest: (i) when the incoming particles 
are a Fermi sea, typically representing the host metal and 
(ii) when the incoming particles are a Fermi sea and an 
excited quasi-particle. The former allows for us to cal- 
culate thermodynamical properties from scattering while 
the latter allows us to compute, in principle, single par- 
ticle S and T matrices. These boundary conditions are 
depicted in Figure [TJ 

Time Dependent and Time Independent For- 
malisms: T = 0: There are two different descriptions 
for scattering processes. In the time-dependent descrip- 
tion, the interaction between the conduction electrons 
and quantum impurity is turned on in the far past, at 
t — to and then adiabatically time-evolved using the 
Hamiltonian 

H = Ho + e^'H;ntO{t - to), (2) 

with Ho describing the free electron bath and Hint the in- 
teractions between the quantum impurity and the incom- 
ing electrons, i.e., between the dot and the leads. Before 
t — to the quantum impurity is decoupled from the elec- 
tron bath and the system is described by an asymptotic 
boundary condition, an eigenstate of Hq which we de- 
note |$o) (at r > the system is described by some den- 
sity matrix po describing decoupled leads and the dot). 




E p 



Filled Femii Sea 



Fermi Sea + 
quasi-particle 

FIG. 1: The incoming particles in the chiral picture corre- 
spond to the particles on the left of the impurity and the 
outgoing scattered particles are on the right. We take the 
convention that the impurity is at a; = 0. Using this conven- 
tion, incoming particles are located on the negative x-axis, 
a; < 0, and outgoing particles to the positive x-axis, x > 0. 
We typically consider two types of open boundary conditions 
for incoming particles. In the first, the incoming particles de- 
scribe a filled Fermi sea. This is useful for calculating thermo- 
dynamics. In the second, the incoming particles are a Fermi 
sea and an excited quasi-particle. This allows us to calculate 
the single-particle 5* and T matrices. 

At later times, as the baths and the quantum-impurity 
evolve, the interaction is turned on adiabatically^^ from 
the state |$o) under the action of the time evolution op- 
erator 

\^it)) = U{t,to)\'^>o) 

U{t,to) - T{cxp{-iJ^ dt'H{t'))} (3) 

We now wish to establish that a steady-state ensues af- 
ter sufficiently long time- long enough that all transients 
die out. For this purpose one must show that the limit 
to —oo exists, free of infra-red divergences. This has 
been shown to be the case for the Kondo modeU^ under 
the condition that the infinite volume limit is taken first, 
i.e. the system is open and the impurity is coupled to 
good thermal baths. The "openness" of the system pro- 
vided the dissipation mechanism necessary for the steady 
state, allowing the high-energy electrons to relax and es- 
cape to infinity. The adiabatic limit 77 — > is taken last, 
allowing the smearing of the bath levels (level separation 
6 ^ l/L to take place turning the poles in the Green's 
f miction into a cut). 

Under these circumstances |^(t)) become time- 
independent and describes a time independent eigen- 
state which we denote Thus, we can also describe 
our state in a time-dependent picture. In the time- 
independent picture time is traded for space and - for 
the chiral unfolded picture- the far past corresponds to 
incoming particles located at distances x ^ and the 
far future to outgoing particles located at x ^ 0. (see 
Figure [T|) . Under both equilibrium and non-equilibrium 
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conditions (e.g. coupling to baths at different chemical 
potentials), the expectation value of any operator O 



(6} = 



WOW. 



(4) 



is time independent. 

A stronger conclusion can be deduced which is central 
to our construction: the state |5'(0)) = Y^)s becomes, 
by the Gellman-Low theoremi^, an eigenstate of the full 
Hamiltonian H, specified by the initial condition |$o) 
that describes the electrons in the far past (x ^ 0). In 
other words , the state Y^)s is a scattering eigenstate 
of the full hamiltonian H = Hq + Hint, satisfying the 
Lippman-Schwinger equation. 
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E- Ho±iT] 



(5) 



with |$o} - the incoming state playing the role of a bound- 
ary condition imposed asymptotically. The scattering 
state can be viewed as consisting of incoming par- 
ticles (commonly taken to be a bath of free electrons) 
described by |$o), and of scattered outgoing particles 
given by the second term in the above equation. Once 
again two elements are required to fully determine the 
system: a hamiltonian, H, and a boundary condition, 
$o), which describes the incoming scattering state far 
from the impurity. Note that previously, in the time- 
dependent picture, |$o) played the role of an initial con- 
dition rather than a boundary condition. As the impu- 
rity is short ranged the scattering state \'^)s must reduce 
to the eigenstate |$o) when all the particles are far to 
the left of the impurity. This gives a prescription for 
calculating scattering eigenstates for quantum-impurity 
problems. We must construct an eigenstate of the full 
Hamiltonian H with the requirement that when all the 
electrons are to the left of the impurity the eigenstate 
reduces to a prescribed eigenstate of Hq describing the 
free decoupled two baths and the impurity. It is worth 
emphasizing that we never explicitly solve ([5]). Instead, 
we directly construct eigenstates of the fuU-Hamiltonian 
with |^')s that are of the form described above. 

Scattering formalism at finite Temperatures: 
The above discussion can be generalized to finite temper- 
atures. Once again there are two equivalent frameworks 
for quantum-impurity problems, a time-dependent and 
time-independent. In the former, we proceed as in the 
zero temperature case. We consider the quantum impu- 
rity and the baths to be decoupled in the far past, at 
t — to — oo, adiabatically turning on the coupling. 
The Hamiltonian is again given by ^ . The change from 
zero temperature is that the system is no longer described 
by a single eigenstate but must instead by described by 
a density matrix. At t — — oo, the quantum impurity 
is decoupled from the electron bath and the system is 
described by the density matrix 



At later times, the system is described by time evolving 
the density matrix po with the time evolution operator 



p{t)^U\t,-^)poU{t,-^) 



(7) 



with U{t, — oo) being understood as the limit U{t, to — > 
— oo). The expectation value of an operator O can be 
calculated in the usual manner by 



(0) = 



Tr(p(t)0) 
Trp(t) 



(8) 



Again a time independent description can be given. 
Now the boundary-conditions for our evolved density ma- 
trix ps is provided by p^: to the left of the impurity, we 
know that the scattering density matrix pa must reduce 
to pq. Further, the finite temperature analogue of our 
zero temperature condition that our scattering state \^) 
be an eigenstate of H is requirement that the density ma- 
trix Pa commute with the full Hamiltonian in the limit 
?; ^ 0. 

Thus, at T > we consider the incoming states, 
{|0, m)}, the complete set of eigenstates of Ho with en- 
ergies Erm distributed with the probability of each state 
given by the Boltzman weight, e"''^™, 

p^^e-P''" ^Y.^-^''"^\<t>,m){<t>M, (9) 

m 

Using ([7]), the time-independent density matrix pa is 

Pa = C/(0,-(X3)p„C/t(0,-oo) 

= ^ e-''^'"i7(0, oo)|0, m){(t>, m\U\Q, -oo) 

m 

= ^e-'3^'"|^,m)(«',m| (10) 

m 

where we have used ([9]) and in the second line we have de- 
fined the scattering state |^',to) = C/(0, — oo)|<I>,m) with 
incoming particles describe by |$,to). The steady state 
physics is captured by the operator ps which describes 
an ensemble of scattering states weighted by the Boltz- 
man factors determined by the energy of the incoming 
electrons. This form for pa is consistent with the require- 
ment that Pa commute with the Hamiltonian and reduce 
to po to the left of the impurity. We can calculate the 
expectation value of an operator as in ([8]) 



Tr(p,0) 
Trps 



(11) 



po = exp {-fiHo) 



(3 = T-\ 



(6) 



III. THE SCATTERING BETHE-ANSATZ 

We have shown in the previous section that the ther- 
modynamic properties can be obtained from scattering 
eigenstates defined directly on the infinite line with in- 
coming boundary conditions imposed by the lead. In gen- 
eral, constructing such eigenstates is a formidable task 
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due to the strong correlations between particles, and is 
only carried out approximately. But for a special class of 
models, many of which have important direct physical ap- 
plication, the many-particle eigenstates can be explicitly 
constructed using the Bethe-Ansatz wavefunction form. 

The Bethe-Ansatz approach has a long history stretch- 
ing back to Bethe's study of the Heisenberg model^S. 
The approach has been typically implemented on systems 
defined with periodic boundary conditions with respect 
to some finite length L. Subsequently the thermody- 
namic limit is achieved sending L to infinity maintaining 
a constant density. If a field theory limit is to be taken, 
then a further scaling (or universality) limit is required. 
By means of this "Traditional Bethe Ansatz" (TBA) ap- 
proach the thermodynamics of several impurity models 
was studied in great detailjiii^ii^. 

Scattering, on the other hand, must by definition take 
place in infinite systems with no boundaries - open sys- 
tems in our terminology. To compute scattering proper- 
ties a different formulation is required. This can be seen 
from several points of view. To begin with, particles must 
come in from asymptotic regions and after scattering oc- 
curs, escape again. Thus the system must be open to 
allow the flow of particles and energy in and out of the 
system. Furthermore, there must be a way to distinguish 
between the incoming particles, typically bare particles, 
eigenstates of Ho but not of H , and the renormalized 
quasiparticles that are the eigenstates of the latter but 
not of the former. Expressing the bare particles in terms 
of the renormalized quasiparticles and vice versa lies at 
the heart of the scattering theory. 

Thus the traditional (or closed) Bethe-Ansatz (TBA) 
and the Scattering (or open) Bethe-Ansatz (SBA) nat- 
urally address different sets of questions. The natural 
questions that can be addressed using the first are ther- 
modynamical. The TBA, by using a periodic system and 
requiring wave-functions to be self-consistent, reproduces 
the full renormalized excitation spectrum of a quantum- 
impurity model. With the knowledge of the spectrum, 
one can use statistical mechanics arguments to calculate 
the thermodynamic quantities. Boundary conditions (pe- 
riodic or otherwise) imposed on a finite length system 
are essential to this approach. But all knowledge of the 
bare theory is lost. As such, the TBA is unable to tackle 
questions about the scattering properties of the quantum 
impurities. Scattering relies on working in systems with 
bare particles and open boundary conditions- namely sys- 
tems of infinite extent with no boundaries. 

There is a price to pay for working in open systems. 
The wavefunctions are no longer normalizable and one 
does not have recourse to thermodynamic concepts such 
as free energy. Thus, while the SBA is essential for ana- 
lyzing scattering properties of quantum impurity models, 
it is more difficult to extract the thermodynamics using 
it. Table [T] summarizes the comparison between the two 
approaches. 





SBA 


TBA 


System 


Infinite 


Finite 


Boundary condition 


asymptotic (open) 


periodic 


Wavefunctions 


used explicitly 


not used 


Thermodynamics 


difficult 


easy 


Scattering Properties 


possible 


not possible 


Nonequilibrium Generalization 


Yes 


No 



TABLE I: Summary of differences between the Scattering 
Bethe-Ansatz (SBA) and Algebraic Bethe Ansatz (ABA). 



A. The Bethe-Ansatz Wavefunction 

The central objective of the SBA is to construct on the 
infinite line eigenstates of the Hamiltonian 



H — Hq + Hint — —i 



dxijl{x)dMx) + Hint (12) 



with the condition that far away from the quantum- 
impurity the incoming sector of the eigenstate reduces to 
a prescribed eigenstate, |<i>), of the free-electron Hamilto- 
nian Hq. As such, any scattering state must have a well- 
defined sense of incoming and outgoing particles, with 
the incoming electrons being to the left of the impurity 
{x < 0) and the outgoing scattered electrons those to its 
right. The state |$) can be any eigenstate of Hq. We 
focus in this paper on the case where |4>) is a Fermi- 
sea of incoming particles. However, many other choices 
are possible. In particular, to calculate S-matrices and 
T-matrices of quantum-impurity Hamiltonian one can 
choose 1$) to be a Fermi-sea with one incoming parti- 
cle above the Fermi-sea (see Figured]). 

The choice of |$) describing incoming particles imposes 
an asymptotic boundary condition on the full scattering 
eigenstate. In general, imposing boundary conditions on 
our scattering states is quite difficult. However, when 
the incoming particles are a free Fermi sea, imposing the 
boundary-condition simplifies greatly. The key to this 
simplification is the observation that the natural basis 
for Bethe-Ansatz wavefunctions is not the Fock basis, 
but a new "Bethe basis" described extensively below. 

The SBA constructs eigenstates of the Hamiltonian us- 
ing wave-functions of Bethe-Ansatz typeSS.. The Bethe- 
Ansatz utilizes the integrability of the Hamiltonian H to 
divide multi-particle scattering events into two-particle 
scattering events characterized by the two-particle S- 
matrices, S'^ derived from H. The integrability of 
the Hamiltonian translates in this language into a self- 
consistency condition on the two-particle S-matrices 
known as the Yang-Baxter Equation (YBE)^'^ ensuring 
that all multi-particle interactions can be consistently 
broken-up into pair-wise interactions. The consistent 
wavefunctions of the Bethe form, which we collectively 
refer to as Bethe-Ansatz wavefunctions, are eigenstates 
of the Hamiltonian. 
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We restrict our analysis to quantum-impurities coupled 
to non-interacting electrons. We further assume that par- 
ticle number is conserved, the Bethe-Ansatz wavefunc- 
tions all have a definite number of particles, N, and there 
arc only local interactions: two particles can interact only 
if they are at the same point. To write a Bethe-Ansatz 
wavefunction, it is necessary to divide the configuration 
space into A^! regions according to the ordering of the par- 
ticles on the infinite line. For example, we can consider a 
region where particle 5 is to the left of particle 7 which is 
to the left particle 9 etc., {x^ < Xj < Xg . . .). Each such 
region is labelled by a permutation Q in the symmetric 
group, Sn+i. Since a particle i and j can only interact 
when they occupy the same position Xi — Xj, there are 
no interactions in the interior of these regions. Within 
each region, the Hamiltonian H reduces to Hq and the 
eigenfunctions are sums of plane waves. The most gen- 
eral wave-function of the above form is (with xq = the 
position of the impurity) 

\BA,{p}) = J dxi.-.dxNe'^^P^''^ (13) 

N 



A = in the region Q = I and the S-matrices , 
\BA,{p}) ^ j dxi.-.dxNe^^P'"'^ (15) 

N 

Y.^SQAU...^J(xQ)\{^i^{x,)\Q). 
Q J=i 

The energy of a Bethe-Ansatz wavefunction (|15p is 
given by -E = J2jPj- Note, however, that the Bethe- 
Ansatz wavefunction with Bethe-Ansatz momenta, {pj} 
is degenerate in energy with all other Bethe-Ansatz wave- 
functions {p'j } with p'- = E = V . Pj . Thus, there are 
an infinite number of degenerate Bethe-Ansatz wavefunc- 
tions of the same energy for any Hamiltonian. Generi- 
cally, a scattering state with energy E' is a sum of Bethe- 
Ansatz wavefunctions (HH) 

l*>= E Qpil^AM), (16) 

with C{p} the amplitude in the scattering state of the 
Bethe-Ansatz wavefunction \BA, {p}) and the sum run- 
ning over all sets of Bethe-Ansatz momenta {p} with en- 
ergy E. 



where ^(xq) = 6l(a;Q(i) < Xq,^2) ■■■Xq(n) < 2^Q(o)) and 
Q runs over all A -I- 1! permutations. The state |0, ao) 
denotes the drained Fermi sea {ipa-ixj)\0) = 0) and the 
state of the impurity. 

When a boundary between two regions is crossed, two 
particles interact (multi-particle interactions forbidden 
by Fermi statistics) and hence the amplitude in the re- 
gions across the boundary are related by a two par- 
ticle S-matrix determined by solving the two-particle 
Schrodinger Equation for the relevant Hamiltonian. The 
amplitude in a region Q, Aai...Q!iv is related to the 
amplitude in an adjacent region, Q', differing from it by 
the exchange of neighboring particles i and j, via the 
S-matrix 5"*^ 

where in the second equality we have used the fact the 
two-particle S-matrix S'*-' acts non-trivially only on the 
sectors of the Hilbert space corresponding to particles i 
and j. In general, the matrix relating the amplitude in 
the region Q = I, defined by (xi < X2 < ■ ■ ■ < xn < xq) 
is related to the amplitude in region Q, (a;Q(i) < xq(^2) < 
... < xq(n) < a;Q(o)), by an S-matrix given by a 
product of two-particle exchange S-matrices 5"^ along 
the path leading from J to Q. Since many paths can lead 
from I to Q consistency requires that 5*^ be uniquely de- 
fined in a path independent way. This is assured by the 
Yang-Baxter condition^^. Thus, the Bethe-Ansatz wave- 
function can be written in terms of a single amplitude 



B. The Bethe-Ansatz Basis 

To construct scattering eigenstates for integrable quan- 
tum models the Bethe-Ansatz wavefunction exploits the 
large degeneracy of the linearized free electron gas. As 
taught in standard chapters on degenerate perturbation 
theory the correct basis of states in a degenerate sub- 
space to perturb from is the one that diagonalizes the 
perturbation, or equivalently, the one to which the sys- 
tem returns once the perturbation is turned off. This is 
precisely the intuition behind the "Bethe basis" of a non 
interacting field theory. A Bethe basis for a free electron 
gas is the basis inherited from the interacting quantum- 
impurity theory when the impurity is removed, or when 
the system is studied far from the short range impurity. 
The basis is defined by the presence of a non-trivial two 
particle S-matrix S^^ between the right moving free elec- 
trons in Hq. Indeed, a moment's reflection shows that as 
the particles move with the same velocity (to the right 
with vp = 1) an S-matrix does not indicate interaction 
but a choice of basis. 

We now discuss the Bethe basis in more detail. For 
a quantum-impurity model, there are two kinds of two- 
particle S-matrices: those that describe electron-electron 
scattering, which we denote 5'^ , and those that describe 
impurity-electron scattering which we denote S'^^ . The 
S-matrices and S'^^ are determined by the impurity 
interaction term Hjat in (fT2|) and the Yang-Baxter con- 
sistency conditions. 

Imagine turning off the coupling to the impurity in 
dni) so that Hint = 0. Then reduces to the free-field 
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Hamiltonian Hq and the electron-impurity S-matrix, S'"-' 
reduces to the identity, S*"^ 1. The electron-electron 
S-matrix 5*-' , however, does not change. This leads to the 
somewhat surprising conclusion that Bethe-Ansatz wave- 
functions of the form with S^^ ^ 1 are eigenstates 
of the free field Hamiltonian H^. 

This can be understood as follows. Consider the first 
quantized version of Hq. In the two-particle sector, the 
first quantized Hq is given by Hq = —i{dxi +9x2)- Notice 
that any wavefunction of the form 



|2;pi,]32;g) 



dxA'' gi{pi+q)xi + {p2+q)x2 



<(a:i)V4(^2)|0) 



(17) 



is an eigenfunction of Hq with energy E — pi + p2 (the 
{ai} label the internal degrees of freedom of the free elec- 
trons). Since q can take on any value, there is an infinite 
number of such states. Any sum of eigenfunctions of the 
above form is also an eigenfunction of the Hq with energy 
E, 



|2;pi,P2) = 



(18) 



J2 yda;idx2e*^'^^^+*P^^M«^„^e^«(-i--^)^t^(a;i)V4(x2)|0) 



where to go from the first line to the second line we have 
used the fact that A%^^^e^'^'^^^~^^^ is the general ex- 
pression for the Fourier transform of an arbitrary func- 
tion, f{xi — X2), of xi — X2- Thus, due to the large 
symmetry of the free electron problem, there is an in- 
finite number of degenerate two-particle eigenstates for 
Hq. The above argument easily generalizes to more than 
two particles: any function of the form 

\N)= J dxi . . .dxNe^"=^P^''^Y[fc,a,{xi - Xj) 

nV4,(^.)|0) (19) 



i<j 



is an A^-particle eigenstate of Hq. Since 0{xq) = 
Y[i<j ^(^Q(i)-Q(j))i of that form [TO] we conclude that 
the most general A^-particlc Bethe-Ansatz wavefunction 
with S'^ a product of electron-electron S-matrices, S^^ , 

\N,BA)^ J dxe'^^P^^^{S'^AU,„^,e{xQ) 

(20) 

is an eigenstates of Hq. However, for 5''^ 7^ 1 (which 
implies S^^ 7^ 1), it is clearly not of the usual Fock-basis 
form, 

\N,F)^ I dfe^^.''^-M„,...„„<(xi)...V'L(^A^)|0)- 

(21) 



The different choices for S*'^ , and in turn S'^ , correspond 
to different 'Bethe-Ansatz' bases for free electrons. The 
choice of S^^ imposed by the impurity interaction corre- 
sponds to working in a particular "Bethe-Ansatz" basis 
for the problem. The usual Fock basis corresponds to the 
choice S"-' = 1. 

We now proceed to discuss the relationship between 
the Bethe basis, with S^^ ^ 1 and the Fock basis S**^ = 1. 
We denote, for a particular choice of a consistent set of 
matrices 5*'^, the resulting Bethe-Ansatz wavefunctions 
by {\BA)n} where n enumerates all possible choices for 
the {p} and Aa^,,,af, in (^0]) . The set of Bethe-Ansatz 
wavefunctions {|-B^)„} form a complete basis for our 
Hilbert space of Hq in the limit of infinite size and par- 
ticle number. In quantum mechanics, different basis for 
the Hilbert space are related by unitary transformation. 
Thus, we can formally define an operator U that relates 
the Fock basis {\F)m} to the BA basis U maps 

states in the Fock basis ([21]) to states in the Bethe-Ansatz 
basis ([^(7| . In general, the matrix U relating the two basis 
is quite complicated since a single state in the Fock basis 
\F)i maps onto a sum of wavefunctions of the Bethe- 
Ansatz form \F)n — > '^^Unm\BA)m- For example, in 
([T9|) we saw that the two particle eigenstate is actually a 
sum over many Fock states of type PT|) . 

However, U simplifies greatly if we restrict ourselves to 
asking how the ground state of Hq in the Bethe-Ansatz 
and Fock basis are related. For a systems with unique 
ground-states, U must map the Fock basis ground state, 
\N,F)gs to the ground state in the Bethe-Ansatz basis 
\N, BA)gs . Thus, the ground state in the Fock basis maps 
to a single wavefunction of the Bethe-Ansatz form (|20p . 
Since the ground state of Hq is a free Fermi-sea, it follows 
that a Fermi-sea can be represented by a single Bethe- 
Ansatz wavefunction. In the sections that follow, we will 
restrict ourselves to this case where we represent a free 
Fermi sea, the ground state of Ho in both basis. 



C. Imposing Asymptotic Boundary-Conditions 

The goal of the SBA is to construct eigenstates of the 
Hamiltonian (jl2p satisfying the asymptotic boundary- 
condition that the incoming particles are a prescribed 
eigenstate, |$), of Hq. We focus on the simplest case 
when incoming particles come from a bath and are a free 
Fermi-sea. Central to the imposition of any boundary- 
condition on the fully interacting Bethe-Ansatz wave- 
functions is the observation that these wave functions 
pick a particular Bethe-Ansatz basis for the free Hamil- 
tonian Hq. Thus, the boundary condition, typically 
formulated in the Fock basis, must be reformulated in 
the natural basis for the scattering state wavefunctions, 
the Bethe-Ansatz basis. The antagonism between the 
Fock basis, natural for boundary-conditions, and the 
Bethe-Ansatz basis, natural for wavefunction is at the 
heart of many of the SBA. We discuss only the zero- 
temperature case. The generalization to finite tempera- 
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tures is straightforward. 

Recall that the incoming electrons in our chiral pic- 
ture are electrons to the left of the impurity, x < (see 
Figure ??). Thus, the asymptotic boundary condition re- 
quires that the scattering state reduce to the eigenstate 
of Ho, 1*} = \^)i,aths (X) \ad), when all particles 

are to the left of the impurity, {xj} < 0, with \^)bath 
a state describing a Fermi sea of free electrons. In gen- 
eral, the scattering state I^P) is a sum of wavefunctions of 
the Bethe- Ansatz form ([T5| . The amplitudes of the dif- 
ferent Bethe- Ansatz wavefunctions C^p. j are determined 
by the asymptotic boundary condition. It was argued in 
the last section that the \^)baths can be written using a 
single Bethe- Ansatz wavefunction of the form (PU)) . Thus, 
in the case where the incoming particles are described by 
\^)baths, our scattering state j^I') can also be described by 
a single Bethe- Ansatz wavefunction. The incoming elec- 
tron corresponds to the regions in the wavefunctions of 
the form 6{xq';xo) = 0{xq>(^^ < xq,(^2) < ■ ■ ■ < xq'(^n) < 
xq) with Q' a permutation of the electrons in the 
problem. Since there are no electron-impurity scattering 
events in these regions, 

SQ' 

can be written entirely in 
terms of the electron-electron scattering matrix S^^ and 
the scattering state j^f) reduces to l^*") when all elec- 
trons are to the left of the impurity. 



1*) 



= / dxi . . .dxN e'^^''^''' 



(22) 



I 



N 



ai ...Qjv„ ;ao 



0iXQ, 



;a;o)n<(^^)|o)- 



The right hand side is precisely of the form ([20|) . We 
therefore conclude that j^*) reduces to eigenstate of Ho 
in the Bethe- Ansatz basis when all particles are to the left 
of the impurity. This leads to the observation that when 
the incoming particles are a free Fermi sea, imposing the 
asymptotic boundary conditions corresponds to choosing 
the amplitude Aq^...q,j^ and the Bethe- Ansatz momenta 
{Pj } foT o, single wavefunction of the form \20(l such that 
4'^} describes a Fermi sea. 

As is usual in the Bethe- Ansatz, we do not seek to 
determine the BA momenta {pj} in the thermodynamic 
limit, computing, instead, the distribution function for 
the BA momenta, p{p). For an infinite system, the distri- 
bution p{p) and the amplitude Aq^ . q^ are independent 
of the procedure used to arrive at them^i. This observa- 
tion allows us to find p{p) and using an auxihary 
Algebraic Bethe Ansatz problem for a system of free elec- 
trons on a finite ring of length L' with Hamiltonian Ho 
and two-particle S-matrices, S]P . In the limit L' oo, 
the distribution function for the BA momenta and am- 
plitude in the auxiliary problem will coincide with those 
of the physical system. p(p) and Aa^...a,^ are obtained in 
the auxiliary problem in the usual way by requiring that 
the wavefunction be periodic. In particular, the ampli- 
tude Aq^ . q^ and the BA momenta {pj} must satisfy the 
auxiliary Bethe- Ansatz equations. 



This program is carried out explicitly for the IRLM and 
Kondo models in the appendix. 

To summarize, the imposition of the asymptotic 
boundary condition on the incoming particle greatly sim- 
plifies in the special case where the incoming particles 
are a free Fermi-sea. The scattering state lip) can be de- 
scribed by a single Bethe-Ansatz wavefunction and the 
imposition of the boundary condition is reduced to find- 
ing the amplitude and BA momenta for this BA wave- 
function. These are found by using the TBA to treat 
the auxiliary problem of free electrons on a finite ring 
of length L' in the appropriate Bethe-Ansatz basis. The 
amplitude and the BA momenta of the auxiliary problem 
coincide with those of the scattering state in the limit 
L' oo. 



D. Computing with Scattering States 

Thus far, we have discussed the explicit construction of 
the scattering states j^*) for integrable quantum impurity 
models. We proceed now to compute the expectation 
values of physical quantities in the scattering eigenstates 
using (g]) 



(*|0|*) 



(23) 



The calculations of expectation values are greatly sim- 
plified because we work directly with infinite systems. 
Technically, this is because for strictly infinite systems, 
we can ignore all but one term in the Slater-determinants 
occurring in the above expression. This simplification is 
the mathematical expression of the physics for infinite 
systems: electrons that scatter off the impurity "leave" 
the system and never return to scatter off the impurity 
again. 

Consider first the overlap between two Bethe-Ansatz 
wavefunctions. Since they are given as a sum of 
plane waves in each region Q, the overlap of two such- 
wavefunctions, {BA,{pj}\BA,{kj}), is (suppressing the 
internal index aj for notational brevity) 



E 



J dxd^ye^^.i'^^^^ 



-^^y^^d{xQ)d{y^) 



N N 

A{Q)Amo\ n ^{ys) n i'^^m (24) 

s=l 3 = 1 

The Fermions field give rise to a Slater determinant 
^ (_1)^9"(S) J dxdye'^^^''''''-P'y'^0{xQ)e{yQ) 



Q.Q.S 



N 



N 



A{Q)A{Q)l[6{xsu)-yj){25) 
i=i 

where S is a permutation of the N particles. Integrating 
over y, we have 



9 



{BA,{p,}\BA,{k,}) 



(26) 



Q,S 




Q,s 

Thus, we see that this expression is the norm of plane 
waves integrated over a region 9{xq). As is usual we 
regularize plane waves by first placing the system in a 
box of size L whose size is then taken to infinity at the 
end of the calculation. This allows us to consider the 
simpler problem of plane- waves 



lim 



L/2 



L/2 



da;e'('=^-P^)^^6l(xi < X2 . . . < xn)- (27) 



It is straightforward to show that the leading order con- 
tribution in L to this integral is /N\ which occurs 
only if the two sets of Bcthc-Ansatz momenta are identi- 
cal {fcj} — {pj}. This is the statement that plane waves 
are 'orthogonal' even on a region 9{xq) for an infinite 
system. Thus, for infinite size systems we can ignore all 
terms in (p6)) where the kj ^ ps-iQ) for all j.. This leads 
to great technical simplifications as we only need to keep 
terms in the sum (|26l) where Q — 1. Similar, simplifica- 
tions occur when computing the expectation value of an 
operator O between Bethc-Ansatz wavef unctions. 



IV. SCATTERING APPROACH TO THE 
RESONANT LEVEL MODEL 

In this section, we will apply the scattering framework 
to a quadratic model, the Resonance Level Model (RLM). 
Despite its simplicity there is much interest in this model 
because it describes the strong coupling physics of the 
Kondo model. It will be shown that our results agree 
with other approaches to this model such as Keldysh or 
Landauer which can be carried out completely in this 
quadratic case. In the next section we shall apply our 
approach to a fully interacting model with strong corre- 
lations. 

The Hamiltonian for the RLM, 

Hrl = Ho + HHLint (28) 

Hq ~ J d.xtp'^ {x)dxip{x) 

HRLint = t{il;\0)d + h.c.) + eddU, 

describes a local level onto which electrons can hop on 
and off. The energy of the level (relative to the Fermi 
energy) is controlled by e^, related to the magnetic field 
in the anisotropic Kondo2i. Notice, that we have already 
projected to one dimension and there are only right mov- 
ing chiral electrons. As explained in the last sections, in 



Free Electrons 

FIG. 2: The scattering state describes a quantum- 
impurity where the incoming particles {{xi < 0}) are a free 



Fermi-sea with N particles. 



the chiral picture, the free incoming electrons are located 
to the left of the impurity, a; < 0, and the scattered out- 
going electrons are to the right of the impurity, a; > 0. 
The RLM serves as a good pedagogical introduction to 
the scattering framework for quantum impurity models 
since it is quadratic model and we will not have to resort 
to the full machinery of the scattering Bethe-Ansatz. 



A. RLM at T = 0: Thermodynamical Properties 

Consider first the zero temperature thermodynamics. 
We must construct a 'in' scattering state, \^)s, describing 
incoming electrons from the host metal scattering off the 
impurity. The scattering state \'^)s is an eigenstate of the 
full Hamiltonian (|28p such that when all the particles are 
to the left of the impurity |^)s reduces to an eigenstate 
|$o) of Hq describing a Fermi see (see Figured]). 

The RLM Hamiltonian conserves total particle 
number. Hence, we can work in a sector of the Hilbert 
space with a definite number of particles, N. Beginning 
with = 1, the most general single particle eigenstate 
is of the form 



dxe'P''gp{x)ip''{x) + Cpd^i ] |0) (29) 



|lp).s = 



Applying the Hamiltonian leads to Schrodinger equation 

dxQpix) + Vep5{x) = pgp{x) (30) 
tgp{0) + edCp = pCp. (31) 

Taking the ansatz that gp{x) is of the form gp{x) = 
A9{—x) + B9{x) and inserting this into the above equa- 
tion, one has, using the regularization scheme 5{x)9{x) — 
i(5(a;)2i, that 



B 
A 



l+i 



2(p-£d) 



(32) 



Thus, the most general single particle eigenstate is given 
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by 



tg{0) t{l + e'^p) 



(33) 



where to get the second equation we have used pip and 
the aforementioned regularization scheme. For future ref- 
erence it will be helpful to define the single particle scat- 
tering state creation operator 

alix) = {e{~x) + e'^''e{x))t(;^x) + 6{x)epd^ (34) 



Since the Hamiltonian ([28)1 is quadratic, a TV-particle 
eigenstate is given by a tensor product of single particle 
eigenstates. The most general iV-particle eigenstate is of 
the form 

N ^ r 

l*> = 1[®\^P^) - n / dxe^^^^^^^alixjm (35) 

Notice that we have not yet specified the momenta {pj} 
of the state. Since we wish to construct a scattering 
eigenstate, these momenta must be chosen to satisfy the 
boundary condition that when all particles are to the left 
of the impurity our eigenstate reduces to an eigenstate of 
Ho describing the incoming electrons of the host metal 
or lead at thermal equilibrium (see Figure [5]) . At zero 
temperature, this means that the scattering state must 
reduce to |$o) = \^)baths ® \4>d), when all particles are 
to the left of the impurity. Here \(j)d) describes some 
impurity state and \^)baths describes a free Fermi sea 

mbaths = j dxi... dxMe'^"=^ p^^^^^x,) . . . ijHx^m 

(36) 

under the additional condition that the momenta of the 
particles {pj} be distributed according to the Fermi- 
Dirac distribution function. Since we are interested in 
the limit where the number of particles goes to infinity, 
it is sufficient to specify the distribution of the momenta 
instead of the individual values of the momenta them- 
selves. 

The single particle eigenstate ([551) consists of an in- 
coming particle, J dx6(—x)e^P^il;''' {x)\0), and an outgo- 
ing scattered wave, J dx9{x)e^'-P^~^^''^ '^p^x)\0) . Since the 
multi-particle scattering state (j35p is a tensor product of 
the single particle state, when all particles are to the left 
of the impurity \'ijj)s reduces to 



I*) 



(37) 



fdxi... dxN n 0(-x,)e*5:f=iP.-.^t(^^) . . . ^\xNm. 

If we choose the momenta {pj} to be distributed accord- 
ing to the Fermi-Dirac distribution, ([55)1 reduces to the 
expression for \^)baths- Hence, our scattering state is 



given by ([35]) with the requirement that the momentum 
distribution of the electrons be chosen according to the 
Fermi-Dirac distribution. 

Expectation values We can calculate the expecta- 
tion value of operators for the RL model using We 
are interested in calculating the dot occupation Ud = 
{d'^d). Since the multi-particle eigenstate (|35p is a ten- 
sor product of single particle eigenstates ([55)1 , it is useful 
to prove some identities about single-particle scattering 
states. We regularize our system, as is usual in scatter- 
ing theory, by placing the system in a box of length L. 
The physical system corresponds only to the L = oo limit 
and finite L properties are not well defined. A straight- 
forward calculations yields (without loss of generality set- 
ting A = 1 in ([55)1) 



lk\lp) = (38) 
e^P<5pfe + —. -— (1 - 5j,k) 



LS, 



'pk 



i{p — k) 



and 



(ik\dU\ip) = 



(39) 



Thus the overlap of states with the same momenta is 
of higher order in L than those with different momenta, 
so that plane waves are an orthogonal basis for infinite 
size systems. In the scattering framework which works 
directly with infinite size systems, it is sufficient to con- 
sider overlaps only of single-particle states with the same 
momenta. 

Consider now the dot occupation. To leading order in 
L, one finds, combining ([39]), ([38]), ((SIJ, and ([J, that the 
occupation is given by 



N 



{rid) ^ -r 



i ^ 11 



Pi 12 



1 



2A 



L ^ A2 



(40) 



with A = t'^/2, where to go from the first to the second 
line we have used the explicit forms of Cp and e^^p . Since, 
we are interested in the infinite size limit N, L oo, we 
can replace the sum by an integral over the distribution 
of incoming electrons which is given by the Fermi-Dirac 
distribution function, 9{ef — p) to yield 



(ud) = / dp9{ef -p) 



2A 



{p, - eaY + A2 



(41) 



We compare this result to the one from the traditional 
Bethe-Ansatz, defined with periodic boundary conditions 
on a ring of length L. In the usual Bethe-Ansatz, one puts 
the system on a circle and imposes the self consistency 
condition that at x = equals at Xj = L^. 
This leads to the B.A. equations. For this model where 
the two-particle S-matrices are trivial, the B.A. equations 
yield for the energy 



2nni 



(42) 
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The {rij} are integers corresponding to the energy of a 
free electron and the {Sp.} the shift in the energies due 
to the impurity. For future reference, define the 'impu- 
rity' energy as Eimp = hmi^oo ^ J2j ^Pj = I dpp{p)5{p) 
with p{p) the distribution that describes the free elec- 
trons in the Bethe-Ansatz basis. From the Feynman Hell- 
man theoreroi^, we know that 



dE 
ded 



den 



IE 



2A 



dpe{E} -p) 



2A 



(43) 



in agreement with the expression we computed using the 
scattering state formalism (PT|) . 

It is helpful to define an operator that directly yields 
the impurity energy using scattering states. This is 
done by considering the overlap of the outgoing scattered 
waves with the unscattered Fermi-sea. Define a state 
$+) that describes a bath of outgoing particles (i.e. all 
particles are to the right of the impurity) 

» n N 
'' i s=l 

Then, we can define impurity energy alternatively in 
terms of an impurity energy operator, Eimp, that acts 
on scattering state 



Eimp\^) = lim -r- log 



L^oo L ^ \^($+|$ + ) 



A straightforward calculation shows that the expectation 
value of impurity-energy operator 

(44) 

agrees with the expression derived from traditional meth- 
ods. The virtue of this operator is that it can be gen- 
eralized in a straightforward manner to all integrable 
quantum-impurity models. This object is closely related 
to the many-body T-matrix for the quantum-impurity 
model. 

RLM at Finite Temperatures: Thermodynam- 
ical Properties Consider now the finite temperature 
case. At finite temperature, T > 0, the system is no 
longer described by single scattering cigenstatc. Instead, 
we must consider a density matrix of the form (jlOp com- 
posed of scattering states weighted by the thermal Boltz- 
mann distribution. Label the set of A'^-particle scattering 
states by the energy of the incoming electrons m)}, 
with m labelling all possible sets of energies for the par- 
ticles pi < P2 < ... < Pn- We expect that these scat- 
tering states are a complete basis for the Hilbert space, 
and indeed find that this assumption reproduces known 
thermodynamic results correctly. 



We calculate the finite temperature properties of the 
RLM using pT|) . The dot expectation value is calculated 
using dot occupation operator fid = d^d: 



Trp, " TrE„e-^^SJ.|V,m)(^,TO| 
Em,n e"^^" ('0, nlV', m) {ip, m\nd\'4>, n) 



-f3E° , 



|(?/',n|V',TO')P 



The above expressions simplify when we work in the infi- 
nite physical L limit since we can keep only leading order 
terms in L. Recall, that m,m' and n are shorthand la- 
bels for the ordered set of energies of the N electrons 
Pi < P2 < ■ ■ ■ < Pn- Thus, if m ^ n, there is at least 
one electron in each state with different energy. Further- 
more, notice that the overlaps of a single particle eigen- 
states ([551) given by ((551) leading order in L only if the 
energies of the two single particle eigenstates coincide. 
Hence, we conclude from ((55)) that the leading order in L 
contribution to the overlap of multi-particle eigenstates 
comes from states where all particles have the same en- 
ergy, or in other words, when the two states are identical. 
Thus, for the infinite L limit, we can set m = m' = n in 
the above expressions to get 



(nd) = J2 



e {ilj,n\hd\ilJ,n) 



^P(n) 



{ip,n\hd\iJ,n) 



- N 



2A 



P(n) 



E 



E "''(^) 



(45) 



P3<^{Pi'}r. 



where P(n) is the Boltzman probability for the state la- 
belled by n. We can now use a standard trick of statisti- 
cal mechanics and replace the sum over all configurations 
by an integral over the average occupancy of a level of 
energy p, N(p), which in this case is given by the finite- 
temperature Fermi-Dirac distribution function, f{p,T). 
This yields 



dpfip,T) 



2A 



{pj - ed)2 + A2 



(46) 



Thus, the effect of temperature is then incorporated by 
requiring that the momentum distribution of the incom- 
ing electrons be chosen according to the finite temper- 
ature Fermi-Dirac distribution for free electrons. This 
expression is in agreement with known results. We will 
see that the above argument is quite general and that the 
effect of temperature can be generically incorporated by 
integrating over finite-temperature distribution functions 
instead of their zero-temperature counterparts. 

An almost analogous calculation using the impurity en- 
ergy operator Eimp yields that the finite energy impurity 
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energy is 

{E,^p) = j dpfip,T)Sp. (47) 

The great limitation of the scattering formahsm is that 
though we can calculate the finite temperature energy, 
calculating the free energy is much trickier. A free energy 
operator can also be defined for these models though this 
is much trickier and will not be discussed in this paper—. 



V. SCATTERING APPROACH TO THE 
INTERACTING RESONANCE MODEL (IRLM) 
THERMODYNAMICS 

In this section, we compute the zero temperature ther- 
modynamic properties of the interacting Resonance Level 
Model (IRLM) within the scattering framework. The 
IRLM Hamiltonian, 

HiRLM = Ho + Hi = -i j dxil}\x)dx'il>{x) + Hi 

^ -i j dx%lj\x)dxtp{x)+t{ip\Q)d + h.c.) 
+U ip^ {0)2P{0)dU + eddU 

describes a local level, d^, onto which spinless electrons 
hop on and off. There is an additional Coulomb interac- 
tion between the level and electrons. We consider only 
the case where > 0, where the level is above the Fermi 
energy of the electrons. Unlike the RLM considered ear- 
lier, this model is no longer quadratic and we must use 
the full Scattering Bethe-Ansatz (SBA) technology to 
construct scattering states. 

We construct the scattering states. They satisfy 
the Lippman-Schwingcr equation and specifying the 
boundary condition on the incoming particles, j^), and 
the Hamiltonian (|48p . uniquely determine the corre- 
sponding scattering state \'^). In this section, we restrict 
ourselves to scattering states where the incoming par- 
ticles are a Fermi-sea at zero temperature |$o)- Such 
scattering states are sufficient to describe the zero tem- 
perature thermodynamic properties of the IRLM such as 
the dot occupation and impurity energy. 

In principle, the scattering formalism can also be used 
to describe quasi-particle S and T matrices. We defer 
these topics to future publications as they require treat- 
ing more complicated boundary condition for incoming 
particles that includes quasi-particle excitations above 
the Fermi-sea. 



A. Construction of the scattering state 

The scattering states for the IRLM are constructed 
using the SBA, directly in open systems of infinite size. 



L —>■ 00. The most general N-particle eigenstate is of the 
Bethe-Ansatz form 

r ^ 
\{p})=A / dxe^^^<^ J|a],^(xj)lO) (48) 

with 

$(p.fc) = tan-i ( _JliPz:J±_\ (49) 

and 5p and Cp given in l|32|) and ([33|— . Note that is 
the operator that creates a single-particle eigenstate ([M]) 
in the non-interacting RLM. The states \{p}) are a com- 
plete set of states in terms of which a particular scattering 
state can be constructed as a linear combination of by the 
set ({p}) determined by the boundary conditions. In our 
case the boundary condition requires that the incoming 
particles look like a free Fermi sea. As discussed pre- 
viously, for this boundary condition a single state \{p}) 
with appropriately chosen set {p} suffices to determine 
\^)s- In more detail, when all the particles are to the left 
on the impurity, {xi} < 0, must reduce to an eigen- 
state of Hq, |$o), describing a zero temperature Fermi 
sea and a decoupled impurity. When all the {xj} < 0, 
the operators {aj,_ [xj)} reduce to {e^P^'^'tp^{xj)} and the 
eigenstate \^)s reduces to 

r ^ 

Thus, we must choose the {pj} in such a manner that 
the above expression describes a free Fermi sea. 

Despite its appearance the expression on the right 
hand side is an eigenstate of Hq. This can be seen by 
applying ho = "^X^iLi^a:^ to the wave function. In- 
deed, since all particles are right mover the scattering 
S-matric 5* = 6***^^''^^^ describes the choice of a Bethe 
basis in the infinitely degenerate energy subspace of free 
electrons. Thus, for {xj} < 0, reduces to an eigen- 
state expressed in the Bethe basis characterized by the 
two-particle S-matrix S = 6**^^'^^^^. This Bethe basis 
is the natural basis for our problem since, as discussed 
previously, degenerate perturbation theory demands that 
we choose the basis for the free electron eigenstates by 
"turning off" the perturbation, in this case the coupling 
to the quantum impurity. It is worth emphasizing that 
the momenta {pj} should coincide with the usual Fock 
momenta of quasi-particles only when U — and the 
S'^ = 1. 

The boundary-condition on incoming particles must 
be implemented in the Bethe-Ansatz basis with a non- 
trivial two particle electron S-matrix S = e**^^'"''^-'. As 
discussed previously, the requirement that the incoming 
particles be a Fermi sea translates in this Bethe basis into 
the condition that in |<f>o) the incoming particles be an 
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eigenstate of Hq of the form ([20)1 

il^\xi)...^HxNm (50) 

with the additional condition that the distribution for the 
BA momenta of the incoming particles, p(j)), satisfy a set 
of free Bethe-Ansatz equations for an auxiliary problem 
of free electrons on a ring of length L' with a two par- 
ticle S-matrix, S = exp (i<i>(p, fc)). These equations are 
derived in the appendix (|A9p and are given by, 



P{p) 
K{p,k) 



— - / dkp{k)K{p,k) 



1 

2t: 

1 d<^{p,k) 
2tt dp 



(51) 



U 

n (j, + k-2ed)^ + i^ip- 



kf 



The desired scattering state, is given by ([48]) with 
the additional requirement that the distribution of the BA 
momenta, p(p), solves the Bethe-Ansatz equation above. 
The simplicity of the equation follows from the fact the 
ground states in the Fock basis and in the Bethe basis 
are unique. This is no longer the case for excited excited 
states. It is also worth emphasizing that (|52p correspond 
to a free Hamiltonian Hq and thus differ from the usual 
Bethe-Ansatz equations for the IRLM^^ in that they con- 
tain no impurity contribution. 



B. Zero Temperature Properties 

Having constructed scattering states, we now use them 
to calculate the thermodynamic properties of the IRLM. 
In particular, we will use scattering states to calculate the 
zero-temperature dot occupation {fid) — (d^d) and the 
impurity energy Eimp defined as using the impurity en- 
ergy operator At zero temperature, Eimp plays the 
role of the free-energy for all dot thermodynamic prop- 
erties. We then show that our results agree with those 
derived using traditional Bethe-Ansatz techniques. 

To calculate the impurity dot occupation we use (jH) 
which yields 



{nd) 



(*|dtd|*) 



(52) 



with 1^1') as in (j48|) . As is usual in scattering theory, 
we regularize our calculations by placing the system in a 
box of size L. Since scattering is defined only for open 
systems, the physical system correspond to the infinite L 
limit and finite L properties are not well defined. From 
the definition of it follows that d'' dajj_^{xs)\0) = 



S{xs)epj''\0). Thus, 



{^\dU\^) =Y,i-'^yA^ J dydxSixs) 



X 62 



I T,i.zj{sgn{xi-Xj)-sgn{yi-yj))'S>{pi,pj) 



N 



\a,,{y,,)d^a]{x,)\0) (53) 



A very similar calculation yields 

(54) 



AT 



n {Q\a,,{y,,)a]{x,)\Q) 



To proceed with the calculation we note that from (|34p . 
one has the relations 



{aj{xj),al(xs)} 



X 5{xs - x-j) + ep.epj{xs)5{xj) 
{d,al{xs)} = epj{xs) (55) 

The right hand side of the first equation has two terms: 
the first term proportional to 5{xs — Xj) comes from the 
anti-commutation of the fermionic field ip while the sec- 
ond comes from d. When calculating ([53]) and ([54]) keep- 
ing only the first term is sufficient to get the leading order 
in L in the dot occupation since the first term contains 
only one delta function where as the second contains two. 
In explicitly open systems where L in infinite, it is suffi- 
cient to treat the anti-commutation relation as 

{a,{x,),al{xs)}^e^^P^-P^^ (56) 
x[0{~x,)e'^'-^-'''^^e{xs)]S{xs - X,) 

Then, the norm to leading order in L is given by 
(*|*) = A'^J2 (-1)'^"" 



-X3)(*(pi,Pj)-*(p^(i),p^(j))) 



N 



S=l 

The integral over y is trivial. As explained in the last 
section, the leading order in L contribution to such an 
integral comes when e*(P=~P''(=))^= = 1 or precisely when 
the permutation a — 1. In this case the integral is per- 
formed trivially and yields ('I'I'I') = A^L^ . An analogous 
calculation using ([53]) yields to leading order in L that 
{^d'^d\^) = A^L^-i Ks?- Combining these two 

results yields 

{^dU\^) 



N N 



2A 



L^^A^ + ip- CdY 
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where we have defined the hybridization A = In 
the, infinite L, infinite N limit, we can replace the sum 
by an integral over the distribution of BA momenta for 
the incoming particles, p{p) given by (jA9[) to get 



ifid) = I dpp[p) 



2A 



+ (p - ea) 



(57) 



We can also compute the impurity-energy using the im- 
purity energy operator 



imp 



- log ^ 



(58) 



where |$^) is the eigenstate of the free bath given by 
([59)1 with the additional requirement that all particles be 
to the right of the impurity: 



r ^ 

s=l 



^. sgn{xi-Xj)iS>{pi,pj) 



(59) 



These correspond to outgoing free Fermi-sea of scattered 
electrons. In this case, 



N 

X n e-'P^'y^'0{y,,){Omy,,)a]{x,m 
= J2 i-iy^"" f dydx 

N 

X X{e^P'-P-^'^^-'e'-^9{x,)5{y,(,)-x,) (60) 

s=l 

Once again the integral over y is trivial and the leading 
order in L contribution comes from when the permuta- 
tion cr = 1 This yields = A^{L/2)^e'^"=i^''^ . 
An almost identical calculation to the one used to calcu- 
late gives ($+1$+) = A'^{L/2)^. Combining these 
results and substituting in (|58|) gives 



(61) 



We can once again replace the sum by integrals over p(jj) 
to get 



Eimp ^ I dp p{p)Sp. 



(62) 



These results can be checked with those arrived at us- 
ing the traditional Bethe-Ansatz (TBA)^. The TBA 
results are almost identical to those from the SBA ex- 
cept that the distribution p{p) must be replaced by TBA 



distributions pi{p) that include a contribution from the 
impurity. 



Eimp = J dp pi{p)Sp 
(fid) = / dppi{p) 



2A 



A2 + (p - edY 



(63) 



Since, as pointed out in^^, the distributions for the TBA, 
pi{p), differs from the distribution from the SBA, p{p), 
by a term proportional to where N is the number 
of particles, in the L, N ^ oo limit, the SBA and TBA 
expressions coincide. 



VI. SCATTERING APPROACH TO THE 
KONDO THERMODYNAMICS 

In this section, we discuss how the scattering Bethe- 
Ansatz could be used to calculate interesting physical 
quantities in the Kondo model. Due to the complex- 
ity of the scattering state for the Kondo model, doing 
concrete calculations requires the generalization of many 
mathematical methods described in the context of spin 
chains. In particular, we discuss the tantalizing possi- 
bility that many of the methods of Maillet, Terras, and 
collaborators^^ can be generalized to the Kondo model 
where they may allow exact calculation of as yet inacces- 
sible interesting physical quantities such as the impurity 
T-matrix. The section starts with a brief discussion of 
the scattering state that captures the thermodynamics of 
the Kondo problem. In the next subsection, we discuss 
a possible mapping between the Kondo problem and an 
auxiliary 'abelian' problem similar to the IRLM model. 
Finally, we discuss how to calculate quantities in this 
auxiliary problem and discuss how this formalism may 
be generalized. We concentrate only on the zero temper- 
ature properties of the Kondo model. The generalization 
to finite temperatures will be presented later. 



A. The Scattering State 

The scattering state for the Kondo model is signifi- 
cantly more complicated than that for the IRLM. These 
extra complications stem from the non-abelian nature of 
the electron two-particle S-matrices in the Kondo model, 
Sl^ — P*-' . This is already evident in the appendix where 
we represent the free-Fermi seas in the Kondo Bethe- 
Ansatz basis. We focus on constructing scattering states 
where the incoming particles are a free Fermi-sea at zero 
temperature. Such scattering states, using a conjecture 
discussed below, allow one to recover the zero tempera- 
ture thermodynamics of the Kondo model using scatter- 
ing states. 

It was shown earlier that for scattering states with the 
asymptotic boundary conditions that the incoming parti- 
cles are a Fermi sea, that the scattering state can be 
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described by a single Bethe-Ansatz wavefunction. The 
most general Bethe-Ansatz wavefunctions is of the form 



N 



•' Q i=i 

(64) 



with S'^ the product of two-particle S-matrices in the 
Kondo model, S^^ = P'-' for electron-electron scattering 
and S"*" = for electron impurity scattering with 

P^^ is the permutation matrix that exchanges the spins 
of particles i and j^^. 

The asymptotic boundary conditions that the incom- 
ing particles be a filled Fermi-sea now reduce to choos- 
ing the Bethe-Ansatz momenta {pj} and the amplitude 
^c(i...ajvQo so that when all the particles are to the left of 
the impurity are scattering state reduces to eigenstate of 
Ho in the Kondo Bethe-Ansatz basis describing a filled 
Fermi-sea. This state, \^)baths is extensively discussed in 
the appendix and is described by a wavefunction of the 
form (|A10[) with Afci...b„ given by (|A13[) and BA momenta 
{pj} of the form — j-^ with rij integers running from —N 
to 0. The amplitude is written in terms of solutions to 
(|A12[) . the spin rapidities {A^}. 

When all particles are to the left of the impurity, the 
scattering state ((M)) reduces to 



I*) 



(65) 



N 



Q' j=i 

with Q' a permutation of the N'^ electrons, 0{xqi \ xq) = 
6{xqi{i) < 9qi(2) < ■ ■ ■ < xq,(^m') < Xo) with a;o the po- 
sition of the impurity. Since reaching the regions Q' in- 
volves no electron-impurity scattering, the is a prod- 
uct of the electron-electron scattering matrix S'^^ — P^^ 
only. If we choose the momenta {pj} and amplitude 
^bi...bjv as in the paragraph above, (??) reduces to the de- 
sired eigenstate of (jA10[) . Thus, the imposition of the 
boundary-conditions follows directly from the represen- 
tation of the filled Fermi-sea in the Kondo Bethe-Ansatz 
basis. 

Summarizing, the full scattering state is described by 
(??) with the additional conditions that Af,j ...b„ be of the 
form (|A13p with the {A.^} solutions to (|A12|) whose den- 
sity is given by (|A16|) and BA momenta {pj } of the form 
with rij integers running from —N to 0. Choosing 
the amplitude and BA momenta in this way ensure the 
scattering state |^) reduces to a state describing a filled 
Fermi sea \^) baths for incoming particles. 



B. Can we map the Kondo to an abelian 
quantum-impurity problem? 

Having constructed the scattering state, the next task 
is to compute quantum-impurity properties using this 



state. This task is significantly more difficult than in the 
IRLM since the amplitude Abi...f,„ is written in terms of 
lowering B operators of the quantum-inverse scattering 
method. These operators do non commute but instead 
satisfy a complicated algebra. This makes it difficult to 
manipulate them^^. For this reason, it is quite desir- 
able to explore the intriguing possibility that the Kondo 
problem is in fact equivalent to an auxiliary quantum- 
impurity problem similar to the IRLM. The central differ- 
ence between the Kondo scattering state and the IRLM 
is that the scattering states for the Kondo problem are 
constructed using non-abelian two particle S-matrices 
where as the two-particle S-matrix for the IRLM is an 
abelian phase. We call models with abelian two-particle 
S-matrices, abelian quantum impurity problems. In this 
section, we conjecture that the Kondo problem can in- 
deed be mapped to a very particular 'abelian' quantum- 
impurity problem. This abelian quantum-impurity prob- 
lem correctly reproduces the thermodynamics of the 
Kondo model. We conjecture that arguments similar to 
those given by Maillet et al will show that the abelian- 
ization of the problem extends to all quantities allowing 
an easy computation of the scattering properties. 

The starting point for the conjecture are the Bethe- 
Ansatz equations for the Kondo model. These Bethe- 
Ansatz equations are derived using the TBA by consid- 
ering a quantum impurity on a finite ring of length L and 
imposing periodic boundary conditions. They are given 

by 



nM A^-l+ic/2 
7=1 A-,-l+ic/2 



(66) 



M 

n 

5=1,55^7 



Aa - A., 



A|5 — A^ — ic 



\^A^-l+ic/2 j \K^+ic/2 J 



with the additional information that the energy of the 
Bethe-Ansatz wavefunction is £' = TlijPj- The A are 
known as the spin rapidity and parameterize the M spin- 
down particles. We also need the log of these equations 
which yields 



2^ 1 ^ 
P3 = ^n, + -^[02(A,-l)-7r] (67) 



7 = 1 



M 



N^e2{K^~l) + iV*6i2(A^) = -27r/^ +^0i(A^ - Aa) 



5=1 



with 0n{x) — —2 tan^^ nx/c and rij and Ij integers com- 
ing from the logarithm and are the charge and spin quan- 
tum numbers respectively. The energy of the eigenstate 
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IS given 



E 



E 



E 



27r 



E x^^^ 



E 

7 = 1 



^ 27r 

Et"^ 



1 

7 = 1 

A/ 



M 

E' 

(5=1 



E 



27r 



iV' 



M 

E 

7 = 1 



^(A7) 



The open problem in this conjecture is how to map op- 
erators in the original Kondo problem to this new abelian 
quantum-impurity problems. Such a mapping has been 
worked out for the Heisenburg spin-chain by Terras and 
collaborators^^. Due to the close analogy of the Bethe- 
Ansatz equations for the Heisenburg spin chain equa- 
tions, we expect that a similar mapping of operators can 
be performed for the Kondo model. If such a mapping 
can be fully flushed out, the SBA should lead to exact so- 
(68) lutions for many impurity properties such as the impurity 
T and S-matrices. 



The first two terms are the energy of a free-electron 
gas in the spin-charge decoupled Kondo basis and the 
last term is the shift in the ground state energy due to 
the Kondo impurity. Previously, we have defined this as 
the impurity energy Eimp- Thus, for the Kondo problem 
we can write (suggestively) 



1 

Eimp — Y E 



(69) 



with (5k(A) = -6'2(A^) = 2tan-^ (2A/c). We can also 
define a phase <i>x(A^— A^) = 9i{K^— Kg) and a function 
fc(A) = L»6l2(A - 1) with D = N^/L. Then, the second 
Bethe-Ansatz equation in (j67p can be derived from the 
equation 



M 

n< 

5=1 



by taking the natural logarithm of both sides. This sug- 
gestive notation is illuminating because the above equa- 
tion is of the form of the self-consistency monodromy 
equation in the TBA that leads to the BAE 



(5^ 



j-i 



y) 



with S^' = g^<s>K(.A,~A,) ^-^g^ 0) and 5^° = (^^k{A-, 

Thus, viewing the A's as a function of the k's, we 
see that the BAE for the Kondo problem could be de- 
rived from another abelian quantum impurity problem 
of M electrons with an electron-electron scattering ma- 
trix given by S^^ — e'*^(^J ~^=) and electron-impurity 
scattering matrix given by S^'^ — e^^'^^^^\ 

The SBA can be applied to this auxiliary quantum- 
impurity problem in the abelian formulation. The scat- 
tering states are analogous to those of the IRLM model 
with {$,(5} {^k,Sk}- a straight-forward construc- 
tion and calculation using the SBA for this abelian prob- 
lem yield the correct Kondo thermodynamic properties. 
This opens up the possibility that scattering properties 
of the Kondo model can be alternatively calculated in 
this abelian quantum-impurity model where manipula- 
tions of the scattering states are much easier. The scat- 
tering states constructed in the last section are unwieldily 
because they are defined in terms of complicated algebras 
found in the ABA. 



VII. CONCLUSIONS 

This paper outlines a scattering framework for 
quantum-impurity models. Generally, constructing scat- 
tering states for interacting impurity models is quite dif- 
ficult. However, if the model is integrable, these states 
can be constructed using the Scattering Bethe Ansatz 
developed in this paper. The SBA correctly reproduces 
the zero temperature thermodynamic properties of both 
the Kondo model and the IRLM. In addition, it raises 
the exciting possibility that the Kondo model may be 
equivalent to an abelian quantum-impurity problem. 

The scattering framework also gives us insight into 
how the Bethe-Ansatz works. The impurity physics in 
any Bethe-Ansatz basis, always looks like single-particle 
impurity phase shifts, 6. This suggests that the Bethe- 
Ansatz basis diagonalizes the lead electrons so that the 
impurity T-matrix is a phase shift. The complexity of the 
problem is shifted from the impurity-electron interaction 
to finding an appropriate basis for free electrons. This 
observation is essential when using the SBA to calculate 
nonequilibrium properties of the Kondo model. We feel 
that this new perspective on the Bethe-Ansatz may lead 
to new physical insights and is worth exploring in greater 
detail. 
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APPENDIX A: THE BETHE-ANSATZ BASIS 
FOR KONDO AND IRLM 

In general constructing scattering eigenstates is a 
formidable task. However, for integrable quantum- 
impurity models such as the IRLM and Kondo Models, 
scattering states can be constructed using a generaliza- 
tion of the Algebraic Bethe-Ansatz and quantum Inverse 
scattering methods. Consequently, the natural basis for 
these scattering states is not the Fock basis, but rather a 
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new 'Bethe-Ansatz' basis. Central to constructing scat- 
tering eigenstates, is the requirement that far away for 
the impurity, the incoming electrons look like a free Fermi 
sea. In this section, we show how to represent a Fermi- 
sea in the Bcthe-Ansatz basis appropriate to the IRLM 
and Kondo models. For these models, the impurity forces 
two-particle S-matrices S'^^ to be2ii2^ 



IRLM 



exp i tan ^ 



UjPi -Pj) 
2(Pj + Pj - 2ed) 



(Al) 



In this appendix we show how to represent free-electrons 
in the Bethe-Ansatz basis for each of these models. De- 
note these two-basis the IRLM basis and the Kondo basis 
respectively. 



Taking the log and multiplying by —i one has 



Pj 



1 ^ 

-Y.'^iPj.Ps 



s=l 



2TTnj 
U 



(A6) 



with Tij an integer. Notice that the amplitude A has 
dropped out of the equation implying that it may be 
taken to be any constant. Notice that the 'free' Bethe- 
Ansatz equations (|A6[) for the BA momenta of iJp in the 
IRLM basis can be obtained form the Bethe-Ansatz equa- 
tions for the full IRLM Hamiltonian (including impurity 
interactions)^. 



(A7) 



1. Free Electrons in the IRLM Bethe-Ansatz Basis 

We first focus on the IRLM. Particles in the IRLM 
are spinless and labelled by their B.A. momenta pj. 
The IRLM Bethe-Ansatz basis has a electron-electron S- 
matrices of the form (|Aip . A Bethe-Ansatz wavefunc- 
tion in the IRLM Bethe-Anatz basis is given by (up to 
an overall multiplicative phase independent of the {xj }) 

\N,BA) ^ A I dx'e'2:,K^.g5E.<,^s«(:r.-2;,)<i.(p.,„) 



(A2) 



with sgn{x) the sign function which is equal to ±1 if a; > 
O/x < 0. In writing the above expression, we have used 
the identity that {e{-x) + e^*6'(a;)) = e^*5*e3*'^9"(^) . 
As discussed in the main text, to find the B.A. momenta 
{pj} we consider an auxiliary problem of free electrons 
living on a finite ring of size L' . In the L' ^ oo the 
momenta of the physical and auxiliary problem coincide. 

We restrict ourselves to the zero temperature case and 
when td is greater than the Fermi-level of the electrons. 
This is the case considered in—. To derive the Bethe- 
Ansatz for the BA momenta distribution functions, we 
must equate the wavefunction for the auxiliary problem 
on a circle when a particle j is at Xj — and at Xj — L' . 
This gives rise to a Bethe-Ansatz condition of the form 

S^^ . . . S^-^'S^-^'A = 5*^'^ . . . S"+^e'P^^' (A3) 

which implies that 

{S^^-^ . . . S'^S'^ . . . S"+^) A = e-'P^^'A. (A4) 

Plugging in the explicit form of the two-particle S-matrix 
for the IRLM from (|Aip . this equation gives rise to an 
equation for the BA momenta {pj} of the form (noting 
that we can cancel A from both sides since it is a con- 
stant) 



by setting the impurity contribution proportional to 
equal to zero. 

As is usual, we will not be concerned with solving the 
discrete version of this equation but instead will solve for 
the distribution function, p{p) describing the density of 
solutions to the equations in an interval (p, p + dp) . It 
is worth emphasizing that such distributions make sense 
only in the limit L' — > oo. In this limit, we can replace 
the sum by an integral to get 



Pj 



dkp{k)<^{pj,k) + ^j^. 



(A8) 



In the usual way, an equation for the zero temperature 
density p{p) is obtained by subtracting the equation for 
Pj from that for Pj+i and expanding in the difference 
Ap = pj+i — Pj which yields^ 



1 

1 d<^{p,k) 



Pip) = I dkp{k)K{p,k) 

K{p,k) 



(A9) 



U 



(k - Ed) 



27r dp 



{p + k-2edY + ^{p~k) 



oiPjL' 



p»E"=i*(p.,p.) 



(A5) 



This equations are valid as long as is greater than the 
Fermi energy of the lead electrons. Though we do not 
do it here, we could also find the distribution of the BA 
momenta at finite temperatures by considering the free 
Thermodynamic Bethe Ansatz (TBA) equations for 
corresponding to the free zero temperature BA equation 
(|A9l) for Hq. 

Summarizing, in the IRLM basis, there is a non-trivial 
two particle S-matrix between free electrons of the form 
(jAip . The presence of this matrix corresponds to work- 
ing in a Bethe-Ansatz basis for the IRLM that is distinct 
from the usual Fock basis. In this basis, the eigenstates 
of i/o are of the form (jA2p with the multi-particle S- 
matrices given as products of two particle S-matrices 
of the form (jAip . For a free Fermi-sea at zero tempera- 
ture, the distribution for the BA momenta, pijp)^ is given 
by (jA9p not the Fermi-Dirac distribution functions. 
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netic vacuum in the space by \lu) — Y[ 
Then, the amphtude ^fci...b„ in (|A10|) is given by 
A(Ai . . . AM)fc,...b„ = B{Ai + tc/2) . . . B{Am + ic/2)\Lo) 



2. Free Electrons in the Kondo Bethe-Ansatz Basis configuration with zero spin we restrict ourselves to the 

sector where M = ^. Let the solutions of the BA equa- 
We now concentrate of the wavefunction of a free- tions with M = N/2 for the ground state be given by 
Fermi sea at zero temperature in the Kondo basis. In the {^f }• Define (3^ = A^" + ic and denote the ferromag- 
Kondo basis, free-electrons have a two-particle S-matrix ^ t^at i i \ -nN ( 1 

gij _ pij -^^here P*-' is the permutation matrix acting 
on the spins of electrons i and p^. The Bethe-Ansatz 
wavefunction for the Kondo Bethe-Ansatz basis is 

l^) = /'^^e^^.^^-"^(^^)^;.;;tAi...;.«^(^Q) (aio) 

<(a;i)...Vl„(a:^)|0) 

with an appropriate product of two particle S- 
matrices P'-*. Note that the amplitude in the region 
Q = 1, Abi...bfj and the choice of BA momenta {pj} are 
still unspecified. We will once again have to choose these 
appropriately by considering an auxiliary problem de- 
fined on a circle of length L' . In the limit where L' — > oo, 
the expressions from the auxiliary problem coincide with 
those for the infinite-size open system. Thus, we can use 
the beautiful quantum-inverse scattering technology^. 

Once again the BA equations for the auxiliary problem 
are derived by equating the wavefunction for the auxiliary 
problem on a circle when a particle j is at Xj = and at 
Xj — L' . This gives rise to a Bethe-Ansatz condition of 
the form 



where the usual spin amplitude notation A 
ten as A 



fci...bA 



(A13) 



is writ- 
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(^j)ai'...aN^ftl---&N — 

(gij-'^ gj'^gjN gjj+i"' 



(All) 



by specifying the position of the M down 
spins and the operators B{Aj + ic/2) are defined in the as 
is usual in the quantum- inverse scattering matrix^^. The 
B are best thought of as generalized lowering operators 
(T~ that lower the spin of particle j. The BA momenta 
{pj} in (|A10|) on the other hand are trivially of the from 
Pj — ^"'^ . For the ground state, it runs from some lower 
cut-off K = —^77^ to the Fermi energy at zero. 

Just as in the IRLM case, we can talk about the density 
of {Xj}, o'o(Aj) rather than the individual A themselves. 
This approach is valid in the limit: N,L ^ go with D = 
N/L held fixed. To get the distribution for the density, 
we take the logarithm of both sides of the second equation 
in ((Il2l) to get 



bi...bN 
ai . . .ajY 



i6i...6j\ 



-ipjL' A 



M 



N''e{2A^) 



-2ttI^ 



We must choose Af,j...6„ such that it is eigenvector for the 
equation ZjA = ZjA with eigenvalue Zj = e~*^^^. Note, 
that in general there are many solutions to this equation. 
We will be concerned with a single eigenvector, namely 
the ground-state. 

A general method called the quantum-inverse scatter- 
ing method has been developed to solve this problem. 
Let TV and M denote the total number of particles and 
the total number of spin down particles respectively. Let 
us denote the spin Hilbert space of particle j by Vj . Let 
(g) Vj be the iV-particle spin-space. The Bethe- 
Ansatz equations for the ferromagnetic vacuum are^l 



A^ 



(A14) 



(5=1 



with 9{x) — — 2 tan^^(a;/c) and I.y an integer. Since we 
are interested in the {Xj} for the groundstate, we set 
M = N/2. We then consider do (A) describing the num- 
ber of solutions in an interval (A, A + dA). Standard 
manipulations yield the equation^^ 



(To (A) 

K{A) 



2c 
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TT 


c2 + 4A2 
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TT c2 + A2 ■ 





dA'cro(A')A'(A- A') 
(A15) 



Zj = X{aj,f3i . . .(3m) 
As - A 



M . 

n ^ 



M 

n 



IC 



7=1 

A, 



A. 



This can be easily solved by Fourier transform to yield 
the equation for the density of {Xj} in the groundstate, 



As - A^ 



1 



A 



7 



N 



(A12) 



<Jo{A) 



1 



N 



2c cosh ^ A 



(A16) 



Each set of solutions to the Bethe-Ansatz equations 
{Xj} corresponds to a different eigenstate of iJo- Dif- 
ferent choices of M correspond to eigenstates of spin 
A^/2 — M . Since, We are interested in the ground state 



To summarize, a Free Fermi sea in the Kondo Bethe- 
Ansatz basis is captured by a state of the form (jAlOp 
with Ah. h„ given by (|XT3)) . the {A^} solutions to (|M2)) 
whose density is given by (jA16p . The BA momenta {pj } 
are of the trivial 
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There arc commonly two meanings of 'adiabaticity' in the 
literature. In the first, adiabaticity means that when com- 
puting physical quantities the limit ?? ^ is well defined. 
This is equivalent to saying that one can actually turn on 
the interaction infinitely slowly, with the limit — + taken 
after the L —> oo limit. This is the definition we use here. 
The second, stronger meaning associated with adiabatic 
turning on, is that the interaction is turned on so slowly 
that there are no level crossings in the spectrum and there 
exists a one-to-one mapping between the spectrum Ho and 
H, now the — > limit is taken before the infinite volume 
limit. We do not use adiabaticity in this stronger sense. 



